Theory of high-symmetry tetramer single molecule magnets 
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We present a microscopic theory of single molecule magnets. From our exact single-ion spin matrix 
elements for four arbitrary spins, we study the single-ion anisotropy of equal spins exhibiting Td, 
D2d, or Cav molecular group symmetry. Each group generates site-dependent single-ion anisotropy. 
For weak anisotropy, accurate Hartree expressions for the magnetization, specific heat, electron 
paramagnetic resonance (EPR) absorption and inelastic neutron scattering cross-section are given. 
For Z?2d, azimuthal single-ion anisotropy leads to the observed Nil EPR splittings. 
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Single molecule magnets (SMM's) have been a topic of 
great interest for more than a decade, because of their 
potential uses in quantum computing and/or magnetic 
storage, which are possible due to magnetic quantum 
tunneling (MQT) and entangled states. In fits to a wealth 
of data, the Hamiltonian within an SMM was assumed 
to be the Heisenberg exchange interaction plus weaker 
global (total, or giant) spin anisotropy interactions, with 
a fixed overall global spin quantum number s.jfl,] MQT 
and entanglement were studied in this simple model. 

The simplest SMM's are dimers.0,0] Surprisingly, two 
antiferromagnetic dimers, an Fe2 and a N12, appear to 
have substantial single-ion anisotropy without any ap- 
preciable global anisotropy. 0, 0, Although the most 
common SMM's have ferromagnetic (FM) intramolecu- 
lar interactions and contain n > 8 magnetic ions,|3 a 
number of intermediate-sized FM SMM's with n — 4 and 
rather simple molecular structures were recently studied. 
C04 and &4 have s — 6 ground states with spin 3/2 
ions on the corners of tetrahedrons. 0, Q A number of 
high symmetry s = 4 ground state N14 structures with 
spin 1 ions were reported. 0, fill flij l Inelastic neutron 
scattering (INS) experiments provided stro ng e vidence 
for single-ion anisotropy in C04 and a Ni4.|2l|lfJ Fits to 
electron paramagnetic resonance (EPR) M4 data assum- 
ing a fixed s were also problematic. |l3| 

Yet there is no microscopic model of FM SMM's in 
which the MQT and entanglement issues crucial for quan- 
tum computing can be understood. To analyze the dif- 
ferences between single-ion and global anisotropy in FM 
systems, we found exact expressions for the single-ion 
spin matrix elements of four general spins, and compared 
global and single-ion anisotropics in the Hartree approxi- 
mation for the magnetization, specific heat, EPR and INS 
transitions for equal-spin SMM tetramers with molecu- 
lar group Td, D2d, and C^ v symmetries. Surprisingly, we 
also found that each molecular group symmetry generates 
site-dependent single-ion anisotropy, and that azimuthal 
Did single-ion anisotropy has a continuous symmetry, ob- 
served as splittings in the N14 EPR resonances. pd| 

We assume four equal spins s\ sit on opposite corners of 




FIG. 1: D 2d (left) and C iv (right) ion sites (filled). Circle: 
origin. Arrows: local axial (left), azimuthal (right) vectors. 



an orthorhombic prism with sides (a, a, c) or on a square, 
with the geometric center at the origin, as in Fig. 1. For 
molecular (site point) groups g — Td, Did, and C± v with 



Y] j, r n = 0, the relative spin positions are 



<± = e+(-l)"/ 2 ±e-(-l) 



(n+l)/2 



(1) 

(2) 

where e„ = [1 ± (— l) n ]/2.Q In tetrahedrons witho = 
Td, c/a — 1, approximately as in C04 and Cr4.[^,U In 
squares with g — C± v , c — 0, as in one Mn4 SMM and NcT 
(with equal total angular momentum j = 9/2). |l2U15l . ll6| 
In prisms with g = D id, c/a > 1, approximately as in 
Ni 4 , a Fe4 and a Cu4,[ll| or c/a < 1. x,y,z are the 
molecular (or global) axes of each tetramer SMM. 

The local (or single-ion) azimuthal vectors satisfying 
all g = d v symmetries are = (—7+ x + -f~y)/^/2, as 
in the right panel of Fig. 1, and the common local axial 
vector is £^ = z. For g = Td,Did, we take the group 
symmetry-satisfying local axial vectors to be = r n /ao, 
as in the left panel of Fig. 1, where ao = \/2a 2 + c 2 , and 
we set c = c/ao and a = a/ao = a/(1 — c 2 )/2. We define 
the local azimuthal vectors from (zf) T — (a,a,c), 



— (1 + c) cos \i + (1 — c) sin fi 
(1 — c) cos fi — (1 + c) sin/i 
2a(sin \x + cos fi) 



(3) 



and yf 



For /j, = Q, x\ ■ z\ 



0. A rotation 



by the arbitrary angle /x about z\ leads to Eq. ©. The 
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other local azimuthal vectors are obtained by ir rotations 
of x\ , y\ about x, y, z. The y a utomatically satisfy all of 
the mirror planes of D 2 d- Although /i is a degree 

of freedom in D 2 d, no fi choice satisfies the remaining 
group operations (rotations by ±2ir/3 about the cube 
diagonals) of Td- Hence, Td only has local axial vectors. 

In the local coordinates of groups g — Td, D 2 d, G^ v , 
the most general quadratic single-ion Hamiltonian is 

n *i = -£(-us«-3D 2 

n 

+J e [{S n -x 9 n f~{S n -y 9 n f]), (4) 

which for these equal-spin, high symmetry systems has 
site-independent J a , J e . 7i 9 si is invariant under all al- 
lowed g symmetries. Unequal spin values or ligands or lo- 
cal structural distortions from the ideal molecular group 
symmetry, as commonly occur in real svstems.|l7j break 
these g symmetries, and also lead to Dzyaloshinshii- 
Moriya (DM) interactions, which vanish for precise 
molecular groups Td, D 2 d, or CA,,. [ia.ll^ Such lower sym- 
metry cases will be discussed elsewhere, [jflj 

To make contact with experiment, we rewrite Ti a si in 
the molecular (x,y,z) representation, 

n 

+ £ J n,a^)i S ^' S n,0}) . (5) 

J° 2d (v) = J Q (c 2 -a 2 ) + 3J e a 2 sin(2 M ), (6) 
J^ip) = Je~c cos(2 M ), (7) 

where a,/3 — x,y,z, and the remaining non- vanishing 
couplings are jf 4 " = J a - J e /2, J%£> y = J e (-l) n /2, 
Jn,l%(p) = ln[Jaa 2 + J e (l - a 2 )sin(»], J%**(ii) = 

■#(/*). J n, 2 yz = -7n J n(V>)> J n(v) = (-l) n a(j a C + 

J e [cos(2/i) ±csin(2^)]). In Eq. ®, {A, B} = AB + BA 
and we subtracted an irrelevant constant. To preserve 
the g symmetries of Ti 9 si , H 9 si contains the site-dependent 
interactions Jnapi^)- The continuous (rotational) sym- 
metry of the local azimuthal vectors in D 2 d leads to 
anisotropy energy strengths in TL g si that are periodic func- 
tions of /j,, as in one-dimensional optical phonon bands. 

We assume a molecular Hamiltonian of H = Ho+H 9 si + 
Tip, where the Heisenberg and Zeeman interactions are 

Ho = -JS 2 /2- 1 B ■ S - SJ 9 (Sf 3 + Sf 4 )/2, (8) 

where 8 J 9 — 0, J' — J, and — J for g = Td,D 2 d, 
and C 4 „ respectively, S13 = S\ + S3, S 2 4 = S 2 + 
#4, S = S13 + S 2 4, 7 = g[iB is the gyromagnetic 
ratio (assumed isotropic, with <? ~ 2), and B = 
B (sin 8 cos cf>, sin 8 sin cf>, cos 8) is the magnetic induction 
at an arbitrary direction (8, </>) relative to the molecular 
coordinates. The time correlation functions of the clas- 
sical analog of this C 4t , TCq model were published. 16] 



The phenomenological global anisotropy interactions 
usually studied in SMM's are 

Hp = —JbS 2 — Jd{Sl — Sy), (9) 

containing axial and azimuthal contributions, 
respectively. [l| They are generally defined relative 
to the global spin principal axes, which for equal spin, 
high symmetry systems are the molecular axis vectors. 

To take proper account of B in TIq, we construct 
our SMM eigenstates in the induction representation by 
(x, y, z) — (x, y, z) ■ M so that B = Bz^ where 

(cos 8 cos 4> cos 8 sin <j> — sin 9 \ 
-sin0 cos</> . (10) 
sin 8 cos 4> sin 8 sin (j) cos 8 J 

A subsequent arbitrary rotation about z does not affect 
the eigenstates. We then set Ti — 1 and write 



s 2 w; m S2i ) 


= s (s + i)\r s ^ S2i ), 


(11) 


s 2 3 |V^ S24 > 


= 51 3 (S13 + 1)I^ 24 ), 


(12) 


s 2 2 4 IC™ S24 > 


= s M ( S2i + 1)|V™ 4 >, 


(13) 






(14) 






(15) 


AT 


= v (s — m)(s + m + 1), 


(16) 



where = + iaSy with a = ±. For brevity, we de- 
note v = {s, m, S13, s 24 , {s„}}, and write \v) = \tp s s ^ S2i )- 
From Eqs. {TTJ to lfT3)l. (^'|7Y |^) = E v $5 v i tV , where 

E^o = - Js(s + l)/2 - 75m 

-W 9 [si 3 (si3 + 1) + s 24 (s 24 + l)]/2. (17) 

We then transform 7i 9 si and Tip to the induction repre- 
sentation, and make a standard perturbation expansion 
for small anisotropy energies { Jj} — (J a , Jb, Jd, Je) rel- 
ative to |J| and jB. |j| The matrix elements of Tip are 
then obtained from Eqs. (|14|) and l|15ll . However, the 
matrix elements of Ti 9 si contain more interesting physics. 

By using symbolic manipulation software for the 
Clebsch-Gordan algebra, we find the single-ion spin ma- 
trix elements with general {s n } — (si, s 2 , S3, s 4 ) to be 

+^. s+1 C m s _ 1 A< s ">;~7 1 -; 

^13 1 13 1 24' ft 24 

•'K.. <(7^ l r: s , ), (is) 

(v'\S n ,*\v) = 6 m ,, m+ JA* m 5 s ,, s r { s s "lf< n s s , 

\ i 13' i 13' i 24' ;, 24 

c 7-j (T , rn a {s n },- s— l,n 

^13' 13' 24' ft 24 

-k ., ,/r '"A;- s , V(i9) 

s 13' 13 ,a 24' a 24 / 
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C™ = v 7 * 2 - to 2 , (20) 

Z?f ,m = ctV( s -^to)(s-ctto- 1), (21) 

+Ss[ 3 , Sl3 4a^ n (s 2A ,s' 24 ), (22) 

^24. s 24 e "^i^3'"( S 13' S 13) 

+^s' 13 ,s 13 e ^/3s2?S4' ( s 247 s 24), (23) 



A 



{««}. 



S 13' S 13' S 24' S 24 



a^3'"(si3> s i3) = t(1 + 6, Si3 ,s 24 )^' 13)Si3 



_^+^i? S 13,S24^ , 

In \^ si,s 3 ,s u s 13 ,s 13 -l 

T1 Si,S 3 ,S %A3 + U' 



(24) 



/3^r( S13 ,si 3 ) = 



i-iy 



Vs,S 13 ,S 2 4^s' 13 ,S 13 



fn ^ u si,S3,s"s 13 ,s 13 - 

I /~lSl3 + l> s 24 JT 

+Cr Sl)S3 _ s <5*i 3 , ai3 +l 

„ MS24 4S— Sl 3 

pSi 3 ,S 24 _ '/si3,Sl,S 3 ^ s 4-s 13 ^s 2 4 

SI ' S3 ' S 4s(s + l) 



(^Sl3:S24 



^z : x : y 



n 4 S 24 4^24 

'/«13 ,Sl,S 3 -"- s -\-s 13 ^1 S -|- S13 _X 

4sV4s 2 - 1 



y/z 2 (Az 2 - 1)' 

x(a;+l)-7/(y + l) 
«(« + !) 



(25) 

(26) 
(27) 

(28) 
(29) 



where 7+ is given by Eq. (0). The prefactors m, ^4g m , 
C™, C™ s _ 1 , Dg' m , and D a _'"L 1 are consequences of the 
Wigner-Eckart theorem for a vector operator. The 

challenge was to obtain the coefficients r^ s "h s '™ / and 

a 13> 13' 24' s 24 

^{s n },s,n ^ ^ Their hierarchical structure based upon 

S 13' S 13' S 24' S 24 

the unequal-spin dimer suggests that analogous coeffi- 
cients with n > 4 may be obtainable. Q 

At arbitrary (0,<fi), the first order corrections E^ x = 
{v\H P + Hsi\v) to the eigenstate energies are 



E v \ — 



Jb ' z (P) [ TO 2 - s ( s +l)]- 



2 

[3m 2 - s(s + 1)] 



cos 2 6 



+^>)sin 2 #cos(2</>)), 



— — ■ _ (a _i_ c 2 _ 2 _ 2 

17 _ ?S,Sl 3l S24 "s,Sl 3 ,S24 '/s+l.Sls.S 



(30) 

(31) 
(32) 



2 (E [(^tS CT+1)As 



<T=±1 



E^Si 3 + (l + (T)/2,S24 



)/ 2 



+ (si 3 <-> S24) ), 



(33) 



<r'=±l ^ 
,0 rf r .s 13 + (l + cr)/2,s 24 

~ l ~ Ly°'si,si,<70-'s+ct(1+ct') 



)/2,s 13 ,s 2 
2 



<T=±1 

+ (S13 <"> S24) ), 



(34) 



where the interactions are given by and following Eq. JJJ , 
and V — {s, S13, S24, si} excludes m. Second order cor- 
rections to the energies, including the symmetry-induced 
site-dependent interactions in Eq. (J5J, will be presented 
elsewhere. [20J For Td, these are important, because their 
site-independent single-ion interactions vanish. 

For Td, D 2 di or C^ v symmetries, E v i has a form en- 
tirely analogous to that of the dimer, |4| except for the 
periodic /i-dependence of the D 2 d interactions, which 
broadens the discrete energy levels into bands, akin to 
one-dimensional optical phonon modes. This and the 
strongly different dependencies upon the quantum num- 
bers V of the coefficients of the single-ion and global 
anisotropy interactions can be employed to provide a 
definitive measure of J, 8J 9 , and {Jj}. This allows 
us to determine how well the phenomenological global 
anisotropy parameters simulate the microscopic systems. 
We therefore present the Hartree approximation in the 
induction representation of four measurable quantities. 

The Hartree approximation is accurate at low ksT/\J\ 
and r yB/\J\ not too small, Q where ks is Boltzmann's 
constant. In this approximation, E v = E v _q + E Vt \ 
is given by Eqs. (fTTjl and (J20J), respectively. We de- 
fine Tr» = Y,u Je"Z7 V v( e ) de - For 9 = D 2d , we set 
E v = E v+ + E v ^ sin(2/i + Xv), where E v ± = (E Vjmax ± 



E u ,min)/2. Then 2?g(e) = \[{e - E^ min ){E^ r 



- £ )]-V 2 



is sharply peaked at £„, max , E V}lnin . For g = T d ,C iv , 
V g v {e) = 5{e-E u ) and E v>mxx . > E v > E v , min . The parti- 
tion function in the Hartree approximation is then Z^jp = 
Tr g u e-^, where /3 = l/(k B T), \ Sl - s 3 \ < s 13 < 8l + s 3 , 

\s 2 - S4 1 < S24 < S 2 + S 4 , \si 3 - S24I < S < S13 + S 24 , 

and — s < m < s. Then, the Hartree magnetization 



Mi 1} (B,T) = 7Tr^(me- /3e )/Z^ i) and specific heat 



C 



g,v 



2 p -f3e\ 



k B f3 2 



Z, 



The Hartree EPR absorption ^x^^S, ui) for clock 



(i) 



(35) 



wise (ct = 1) or counterclockwise (u 



T) circularly 
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polarized oscillatory fields normal to B is 
3 x 9 j(i) = St^T^l ie -^\M vy \ 2 

x [S(e - e' + uj) - S(e' - e + uj)] , (36) 

where M u y = A° m 8 m * , m+< A' ..^.i,,^,^ andTr^, = 

Si/' /e"' ^'( e ')<^ e '- The resonant inductions for g = 
7d, C4„ are strong and centered at 

7 S« = ±w + ^^((l-3cos 2 ^f(0) 

-3sin 2 0cos(20)J^ y (V4)). (37) 

For 5 = Z?2d, the weak center of each B^H, Eq. (|3Tjl . 
is equally surrounded by two strong resonances split by 
7AB 1 ^ 1 S ) = 2\aJ e (z v > - z v )\, where z v = ^xl + y*,y v = 
ica F [3m 2 - s(s + 1)] sin 2 9 cos(20), and = |a 2 {26 37 + 

a F [s(s+l) sin 2 6>+m 2 (3 cos 2 0-1)]}. For = 0, TAB r ( ei = 
6a 2 |(2r7i + <r)a-J e \, nicely fitting s = 4 (a- = 1/7) data 
on three N14 compounds. For c/a w 1.1, | J e | w 0.07K, 
0.63K, and 0.91K.[ll) More B rcs (6,(f>) data are urged. 
The Hartree INS cross-section Sg (B, q, u>) is 

i/' 5, _^ n,n ; — 1 

xe^^-""')^!^^!^)^!^^), (38) 

where a,/3 = x,y, z, g 5 = sin06 >g cos^ q , % = 
sin9b t qSiri4>h,q, and qz — cos6i, iq , Q^^ and <fib,q describe 
the relative orientations of B and q,[4| the r„ and 
(^'I'S'n.fil^) are given by Eqs. Q, ifTHl) . and l(T5)l respec- 
tively, and the scalar q ■ (r„ — r n i ) is invariant under the 
rotation, Eq. (JTJJJ. As for the dimer,Q additional EPR 
and INS transitions with amplitudes higher order in the 
{Jj} are obtained in the extended Hartree approxima- 
tion, but will be presented elsewhere for brevity. |20| 

We presented a microscopic theory of high-symmetry 
single molecule magnets. From our exact single-ion spin 
matrix elements for four general spins, we studied the 
most general quadratic single-ion interactions in equal- 
spin tetramers with molecular group symmetry T^, Z?2d, 
or Ci V . Each group introduces site-dependent single- ion 
interactions, and azimuthal anisotropy with D2d symme- 
try has a continuous symmetry observed as splittings in 
N14 EPR resonances, providing a direct measure of the 
microscopic interaction. We used the Hartree approxima- 
tion to provide explicit expressions for the magnetization, 
specific heat, EPR absorption, and INS cross-section, 
valid at low temperatures and sufficiently large magnetic 
fields. Our procedure is extendable to systems with lower 
symmetry, higher-order single-ion interactions, and pos- 
sibly to systems with n > 4 spins. 
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